8DA068848 


MRC  Technical  Summary  Report  1902 

ERROR  ESTIMATES  FOR  THE  APPROXIMATION 
OF  AN  UNKNOWN  CONSTANT  COEFFICIENT  IN 
A PARTIAL  DIFFERENTIAL  EQUATION 

Richard  S.  Falk 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 

610  Walnut  Street 
Madison,  Wisconsin  53706 


December  1978 


(Received  September  11,  1978) 


\ 


Approved  for  public  release 
Distribution  unlimited 


Sponsored  by 

U.S.  Army  Research  Office  National  Science  Foundation 

P.O.  Box  12211  Washington,  D.C.  20550 

Research  Triangle  Park 
North  Carolina  27709 


UNIVERSITY  OF  WISCONSIN  - MADISON 
MATHEMATICS  RESEARCH  CENTER 


ERROR  ESTIMATES  FOR  THE  APPROXIMATION  OF  AN  UNKNOWN 
CONSTANT  COEFFICIENT  IN  A PARTIAL  DIFFERENTIAL  EQUATION 

Richard  S.  Falk 

Technical  Summary  Report  1902 
December  1978 

ABSTRACT 

A priori  error  estimates  are  derived  for  the  approximation  of  an 
inverse  problem  in  which  it  is  desired  to  identify  an  unknown  constant 
coefficient  in  a one  space  dimensional  parabolic  or  two  space  dimensional 
elliptic  partial  differential  equation  whose  general  form  is  known. 


AMS  (MOS)  Subject  Classification  - 65M30 

Key  Words  - Inverse  problem,  identification  problem 

Work  Unit  Number  7 - Numerical  Analysis 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0Q24 
and  the  National  Science  Foundation  under  Grant  No.  MCS 78-02737. 

f 


SIGNIFICANCE  AND  EXPLANATION 


In  this  paper  error  estimates  are  derived  for  the  approximation  of 
an  inverse  problem  in  which  it  is  desired  to  determine  an  unknown  con- 
stant coefficient  in  a one  space  dimensional  parabolic  or  two  space 
dimensional  elliptic  partial  differential  equation  whose  general  form 
is  known. 

This  type  of  problem  has  a wide  variety  of  applications  (e.g.  in 
pollution  problems  in  rivers  and  oil  flow  problems)  in  which  one  is  try- 
ing to  build  a mathematical  model  to  study  the  underlying  phenomena,  but 
is  not  able  to  measure  directly  certain  material  properties  which  are 
necessary  to  forming  the  model  (e.g.  the  diffusivity  in  a heat  flow 
model ) . 


' 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


ERROR  ESTIMATES  FOR  THE  APPROXIMATION  OF  AN  UNKNOWN 
CONSTANT  COEFFICIENT  IN  A PARTIAL  DIFFERENTIAL  EQUATION 

Richard  S.  Falk 

1.  INTRODUCTION. 

In  this  paper  we  wish  to  present  further  results  on  a priori  error  estimates  for 
the  approximation  of  inverse  problems  in  which  it  is  desired  to  identify  an  unknown 
coefficient  in  a differential  equation  whose  general  form  is  known. 

In  two  previous  papers  we  considered  the  two  point  boundary  value  problem: 

(1.1)  (a  ~ ) + c(x)ua  * f(x),  0 < x < 1 

(1.2)  ua(0)  = g ua(l)  = g2  . 

Assuming  f(x),  c(x),  g^,  and  g^  are  known  the  problems  were  to  determine  a constant 
a and  a function  ua(x)  satisfying  (1.1),  (1.2)  and  the  auxiliary  condition 

. a 

(1.3)  " a dJ  (0)  = g3  (see  l3]) 

or 

(1.4)  ua(x*)  = g4  (see  (41), 

where  g^  and  g^  are  given  constants. 

The  purpose  of  this  paper  is  to  show  how  the  techniques  employed  in  |3)  can  bo 
generalized  to  derive  error  estimates  for  the  approximation  of  an  unknown  constant 
coefficient  in  one  space-dimensional  parabolic  and  two  space-dimensional  elliptic 
partial  differential  equations. 

In  particular,  we  will  consider  the  initial  boundary  value  problem 

(1.5)  ua  - au®^  + c(x,t)ua  = f(x,t),  0 < x < 1,  t>0 

(1.6)  ua(0,t)  =■  g (t),  ua(l,t)  = g,(t),  t > 0 

5 6 

(1.7)  ua(x,0)  * g7<x),  0 < x < 1 
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Assuming  f*c,g,g,  and  g are  known,  the  problem  is  to  determine  a constant 
do  7 

a 

a and  a function  u (x,t)  satisfying  (1.5)— (1.7)  and  the  auxiliary  condition: 

(1-8)  -aua(0,t*)  = g (q  and  t*  > 0 given)  . 

X o 8 

We  shall  henceforth  denote  by  Problem  (P)  the  problem  (1.5)— (l.B) . 

We  will  also  investigate  the  elliptic  boundary  value  problem 

(1.9)  - aAua  + c(x)u3  = f(x),  x f SI 

(1.10)  u3  = g on  F 

2 

where  ft  is  a smooth  bounded  domain  in  K with  boundary  T.  Assuming  f,c,  and  g 
are  known,  the  problem  is  to  determine  a constant  a and  a function  u3(x)  satisfying 
(1.9)- (1.10)  and  the  auxiliarly  condition: 

3 a 

(1.11)  a (x*)  = g*  (g*  given)  . 

on 

Problem  (1.9)— (1.11)  will  henceforth  be  denoted  by  Problem  (E) . 

As  in  (3),  our  main  concern  will  be  to  derive  an  a priori  estimate  for  the  error 
| a - a | where  a is  an  approximation  to  a obtained  by  an  appropriate  numerical 
scheme.  To  do  so  we  will  first  need  to  determine  conditions  on  the  data  under  which 
Problems  (P)  and  (E)  are  well  posed.  These  results  are  contained  in  Sections  3 and  5 
respectively. 

We  then  consider  in  Sections  4 and  6 respectively  approximation  schemes  for  the 
two  problems  and  derive  a priori  error  estimates.  Finally  for  completeness,  we  include 
an  appendix  in  which  we  give  proofs  of  several  lemmas  which  we  need  in  the  derivation 
of  our  main  results.  For  the  most  part  they  only  involve  minor  modifications  of 
results  already  appearing  in  the  literature. 

References  to  much  of  the  work  on  the  determination  of  parameters  in  elliptic 
and  parabolic  problems  can  be  found  in  the  bibliographies  of  [10]  and  111).  In  par- 
ticular we  wish  to  mention  the  work  of  J.  R.  Cannon,  since  some  of  the  results 
obtained  and  techniques  used  here  are  extensions  of  [2]. 
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2 NOTATION  AND  PRELIMINARY  RESULTS. 

2 

For  *2  a bounded  domain  in  ® or  I ■ 10,1)  and  k a nonnegative  inteqer  we 
shall  denote  by  W ' the  usual  Sobolev  spaces  of  function  defined  on  i)  ot  I (the 


intended  domain  will  be  clear  from  the  context)  with  norms 


Ml  k (,  f l|DJ»||[  l"1’ 

“ ,p  Mji-0 


! < p 


Ml  k * " I I|d3u||  . 

wk'  | j I *0  L~ 


k k 2 1 

We  further  denote  by  H the  space  W ' and  by 


the  closure  of  in  the  H norm.  We  will  use  the  notation  llu!lk  to  denote  the 

K 

norms  in  H , and  ( • , * ) to  denote  the  inner  product  in  i)  or  I.  Also  for 


X a normed  space  with  norm 


and  u ; |0,T]  ■*  X,  we  define 


II u II  k (x)  " / llu(-.t)  ||x dt  and 

IMIL  (X)  ’ SUp  Hu(‘  't)  Hv  • 

«o  0<t<T 

Finally  we  shall  let  IT  - (0,1)  * 10, T),  C denote  a generic  constant  depending  only 
on  the  data  of  Problem  (P)  or  (E),  and  C(b#,b*)  denote  a generic  constant  which  may 
depend  on  b#  and  b*  as  well  as  the  data  of  Problem  (P)  or  (E) . 

We  now  state  several  lemmas  which  will  be  needed  in  the  proofs  of  our  main  results. 
The  first  two  are  results  from  perturbation  theory  and  we  provide  proofs  in  the  appendix 
for  the  special  cases  we  consider  in  this  paper,  although  results  of  this  general  kind 
can  bo  found  in  the  literature  (see  for  example  (5>)  and  )7)). 


Lemma  1 : Let 


be  the  solution  of 


a a 

az  + cz  « F,  0 < x •.  1,  t > 0 

XX 


za(0,t)  « 0,  za(l,t)  - 0,  t > 0,  ami 


z (x,0)  » G(x)  (with  G (0)  - G(l)  - 0)  , 


and  let  ; bo  the  solution  of 


0 0 

z ♦ cz  - F,  0 < x < 1,  t > 0 


z (x,0)  - G (x)  . 
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Then  if  F,c,  and  G are  sufficiently  smooth  we  have  for  all  0 < a < 1/4  that 


(2.1) 

(2.2) 

and 

(2.3) 


lza  - Z°1  (*,t)  II  0 <_  K^a1/4 


[\  - zt,l*'t)ll0i  v1/4 


{F  - cz3  - z^)  (-,t)  || Q <_  K^a1/4 


where  K^,  K^,  and  K*  are  constants  depending  only  on  the  data  c,F,G  and  t. 

4 5 6 

a 

Lemma  2 : Let  z be  the  solution  of 

a a 

- aAz  + cz  * f - cot  in  0 
za  = 0 on  r . 

Then  if  f,c,ui  and  T are  sufficiently  smooth,  there  exists  0 < < 1 such  that 

for  all  0 < a < e“?,  ||  f - cut  - cza||  < Ca*^8  for  some  constant  C depending  only 
~ 0 l4 
on  f,c,ut,  and  (1. 

The  third  lemma  is  a statement  of  the  maximum  principle  for  parabolic  equations 
which  we  shall  use  in  the  following  form.  (A  general  reference  for  results  of  this 
type  and  for  similar  results  for  elliptic  problems  which  we  shall  also  use  is  [121). 
Lemma  3:  Let  v satisfy 


v - av  + c(x,t)v  = F(x,t) 

t XX 


v(0, t)  = G1(t),  v ( 1 , t ) = G2(t) 
v(x,0)  = G (x) 

where  G^  (0)  * G^tO),  G2(0)  ■ G^d)  and  c.F.G^G^  and  G3  are  smooth  functions. 

If  a > 0,  c ( x , t ) ^ 0,  F (x,  t)  ^ 0,  G^  ( t)  < 0,  G^(t)  ^ 0,  and  G3<x)  <_  0 , for 

(x,t)  t (0,1)  x (0,T) , then  v(x,t)  0,  (x,t)  f 10,1]  x [0,T], 

We  remark  that  even  if  the  compatibility  conditions  G^IO)  = G^JO),  G2<0)  = G^d) 
do  not  hold,  we  still  have  v(x,t)  £ 0,  (x,t)  e (0,1)  x (0,T],  even  though 

v 1 C°[IT).  To  see  this  take  a sequence  G^fx)  <_  0 satisfying  G^tO)  = G2(0), 

G^(l)  • G2(0)  and  converging  to  G^  in  L2(0,1).  Denoting  by  vn(x,t)  the  solution 
of  the  corresponding  initial  boundary  value  problem,  we  have  by  standard  energy 


I jri 


i . J 


-4- 


estimate*  that  v’N’.t)  converge*  to  v(*,t)  in  L ^ (0,1).  Now  vn(x,t)  ^ 0 by 
lemma  3 and  v1'  and  v are  both  continuous  functions  for  t > 0.  Hence  if  v(x,t)  ' 0 
for  some  (x,t)  * (0,1)  * (0,T],  then  we  could  find  a function  i|'  > 0 such  that 
( v ) <_  0 for  all  n,  but  (v,i}>)  > 0.  This  contradicts  the  convergence  of  v"  to 
v in  L.,  (0, 1) . 


* 


3.  WELL  POSEDNESS  OF  PROBLEM  (P) . 


In  this  section  we  wish  to  determine  some  simple  conditions  under  which  problem 

(1.8)  is  well  posed,  i.e.  there  exists  a unique  positive  solution  to  the  equation 

Q(a)  = -au  (0,t*)  = g which  is  continuously  dependent  on  g (we  shall  assume  that 
x o 8 

c,  and  g are  fixed), 

b 6 7 

Although  theorems  under  alternative  hypotheses  are  possible,  we  shall  assume  that 


the  data  c. 

f,g5'q6'  an<^  g7  dre  sufficiently  smooth 

conditions: 

(HO) 

g5(0)  = g?  (0) , g^(0)  = g?(l) 

(HI) 

c(x,  t)  _>  0,  (x,  t)  f I 

(H2) 

g^(t)  :>  0,  t e 10, TI 

(H3) 

g’(t)  <0,  t e [ 0 , T ) 

(H4) 

b — 

g’7  (x)  <_  0,  x e (0,11 

(H5 ) 

f ( x , t ) - c(x,t)g7(x)  <^0,  (x,  t)  e I 

(H6) 

c ( x , t ) - 0,  ( x , t ) e I 

(H7) 

t — T 

f ^ ( x , t ) - ct(x,t)g7<x)  - c (x, t) { (1  - x) [g5 

(H8) 

g*(0)  < o 

(H9) 

t* 

/ g’ (T)dT  - q’  (0)  > 0 . 

We  note 

0 

1 that  these  conditions  imply  that 

(3.1) 

gc(t*)  > q (t*)  . 

5 6 

To  see  this,  write 

g6<0)  = g?  ( 1 ) = g?(0)  + g^(0)  + g"(f,)  = q^  (0)  + g}(0)  + g^(U 

(using  (HO)  and  a Taylor  series  expansion).  Hence  by  (H4)  and  (H8)  , gf(0)  ^_g,.(0) 
with  strict  inequality  holding  if  g^O)  < 0.  By  (H2)  and  (H3)  our  desired  result  is 
true  unless  both  g7<0)  = 0 and  g^.(t)  =0,  t f [0,t*l.  (H9)  rules  out  this 

[>ossibi  1 i ty . 


I 
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We  now  prow  the  following: 

Theorem  1;  Suppose  that  hypotheses  (H0)-(H9)  are  satisfied.  Then  for  any  g > 0, 

B 

there  exists  a unique  positive  solution  of  C(a)  5 -aua(0,t*)  * g . Furthermore,  for 

X 8 


all  0 v a,b  v a* , 
(1.2) 
whe  re 


| a - b | <_  |o(a)  - Q(b)  |/K(a‘) 


K(a*)  «■  - j 


qMr) 

J; dr  - g^(0)  . 


0 rna* (t  - t ) 

hwt : Wnte  u'*  - z3  + r where  r - (1  - x)qc(t)  + xg,  (t)  and  zd  satisfies 

3 6 


(3.  J) 


a a a 

z - az  + cz  » F,  0 < x < 1,  t s 0 
t xx 


za(0,t)  = 0,  zd(l,t)  =0,  t > 0 


z (x,0)  = G(x),  0 < x < 1 

when'  F f - er  - r and  G = g^  - r(x,0). 

The  condition  v'(a)  _ -au'*(0,t*)  = q then  becomes 

x 8 


(3- 4)  £(a)  _ - az'  (0, t*)  + alqc  (t*)  - g (t*)  ] = n . 

X DDO 

By  integrating  (3.  )) , we  get 


azd(x,t)  - az'*(0,t)  = -/  [F  - czd  - za]  (s)ds 
x x I t 


Integral  inn  again  a:ui  applying  the  boundary  conditions,  we  have 

1 x 

az'(0,t)  = / / [F  - cza  - zd](s)ds  dx  . 

* 0 0 1 

Applying  the  Schwarz  inequality  and  setting  t » t*,  we  easily  obtain 

2 


( V S) 


!azd(0,t*)  | < ~ ||  [F  - cz3  - zd)  (•  ,t*)  ||0  . 


It  then  follows  that  l im  = 0 since  1 im  ||  |F  - cz'*  - ?.'*  1 ( - , t ■*  > ||  0 by  lemma  1. 

a Hi  a »0  t 1 

We  will  now  show  using  a standard  energy  argument  that  j|  z'*  ( • , t*)  ||  ^ and 

j z (*,t*)  ► 0 as  a -*  ®.  ?his  coupled  with  the  fact  (3.1)  that  q it*)  - g (t*)  " 0 

t U 3 t > 
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will  prove  that  lim  Q(a)  = ■*■<».  The  existence  of  a positive  solution  to  Q(a)  = g 

a KO  f 

for  positive  g will  then  follow  from  the  continuity  of  Q(a). 

O 

Multiplying  equation  (3.3)  by  za  and  integrating  by  parts,  we  have 


1 d 


? fz  |1  Z” ( - . t ) ||"  + all  z (-,t)  |r  + ||  (►'cz)  ( - , t ) 


dt 


•0  M x ' ' "0 

a . 


<F*Z  > ±l|F(.,t)||0||za(.,t)||0 


i6l  ll*‘(-.t)|lJ^II«aC.t,||J 

(using  the  Schwarz  and  arithmetic-qeometric  mean  inequalities).  Since  za(0,t)  = 0, 
we  easily  get  !j  z^  ( • , t ) ~ 2 ||  z‘l  ( • , t)  ||  ~ and  so 

Hence 


_d_ 

dt 


■at  IU*lst)||?] 


0,_^.“llrc.t>|i; 


and  after  integrating  we  get 


l*a(s tlllji  .-,t||za(-,0)||^i  1 .a(s-t>||F(.  ,s)  ||*  ds  . 


Since  t*  » 0,  we  obtain  lim]|z  (•»t*)||  = 0. 

a 

To  obtain  a similar  result  for  z ‘t , we  differentiate  (3.3)  to  obtain  the  following 
initial  boundary  value  problem  for  z^  (recall  that  G(0)  = G(l)  = 0 by  (HO)). 


( 3.6) 


r ai  r ai  . a a _ 

lztJ  " alztJ  + czt  = "ctz  + Ft 


za(0,t)  = 0,  za(l,t)  = 0 


z^lx.0)  = aG"  - c(x,0)G  + F(x,0) 


Using  the  same  argument  as  in  the  previous  case,  we  obtain 

t 

a (s-t)| 

« 

0 


II  o — e at  !lzt(‘'0)  lip  + i / ea(S"t)||  [-ctza  + Ft]  <-,s)  ||‘ds 


It  easily  follows  that  lim  ||  ( • , t*  ) ||  ^ = 0. 

a ■*** 


-a- 


To  obtain  a continuous  dependence  result  it  is  convenient  to  rewrite 


a a a . . . 

u « y + w + q_,,  where  y satisfies 


yd  - aya  * 0,  0 < x < 1,  t > 0 

t xx 

yd(0,t)  » g^(t)  - (0),  y-'  ( 1 , t ) 0,  t ' 0 


y (x,0)  * 0,  O'  xl , 


and  w'  sat  is  fit 


and  so 


d 3 d 3 „ . 

w - aw  ♦ cw  ■ t - cy  + aa " - eg  _ , 0 «•  x < 1 , t > 0 

t XX 

wa(0,t)  ■ 0,  wd(l,t)  - q (t)  - q (0),  t > 0 

6 6 


w lx,0)  = 0,  0 < x < 1 . 


Q(a)  w -au*(0,t*)  --  -ay‘l(0,t’*)  - .iwa(0,t*)  - uq'(0) 
x xx 


3 a 3 a 

O' la)  ■ -r~  l -ay  (0,t*)l  + — l - aw  (0,t*)]  - q'l 

Ja  x .hi  x 


Now  it  is  not  difficult  to  show  th.it 


-r-  [-awd(0,t*)l  » -a  t—  w"1 10,  t * ) - wa(0,t*l  - [a 

3a  L x J 3a  x x 3x  3a 


.'w  ai 


(0, t *) 


where  satisfies  the  partial  differential  equation 


3 r3w  i ,'w  i .'w  i 

3t  Ua  J " a „ 2 1 3a  ' c Ua  J 


3 r 3w‘ 


\ 

a 3v 

w - c r-1 — > v’, 

XX  vM 


subject  to  homogeneous  initial  and  boundary  data,  and  -----  satisfies  the  partial 


d i t f e rent i a l equat ion 


(3. 10) 


JL  [hi]  _ a J1  [hi] 

3t  L3a  J . 2 l3a  J 


a 1 si 

v = — y 
xx  si  t 


also  subject  to  homoqeneous  initial  and  boundary  d.it.i. 


~ a 

a dw  a 

Setting  p * a + w , we  have  that 


(a)  pa(0#t*)  + ~ l-aya(0#t*)]  - g\(0) 

dX  r ?a  X 


-a- 


+ ag^  + wa 


where  p satisfies 


1 1 1 1 \ a a a a 3y 

(3.11)  p - ap  + cp  = aw  - ac  -r*— 

t xx  xx  3a 


a a 

aw  + cw 

XX 


a . a . Dy 

= w + cw  + aq"  - ac  tt — 
t / i)a 

with  p3 (0 , t ) = 0,  pa(l,t)  = g (t)  - g (0),  and  pa(x,0)  = 0. 

b 6 

We  now  show  that  p(x,t)  <_  0 for  all  0 <_  x <_  1,  t >_  0.  ^ince  pa(0ft)  = 0, 

this  will  imply  that  - — p (0,t*)  _>  0-  To  do  so  wo  observe  that  by  hyjx>theses  (113) 

and  (H4)  g (t)  - g (0)  < 0 and  ag!!  < 0.  Hence  the  desired  result  will  follow 
b b — 7 — 

immediately  from  the  maximum  principle  (lemma  3)  if  we  can  show  that  wa  0,  w‘*  <_  0, 

and  > 0,  (x,t)  e (0,1)  x (0,T) . 

Now  usinq  hypotheses  (H1)-(H5)  it  easily  follows  from  the  maximum  principle  that 


y >_  0 and  hence  that  w*"  <_  0 for  all  (x,t)  € I To  show  that  w^  <_  0 
3y  a 

~ — > 0 we  observe  that  w satisfies  the  initial  boundary  value  problem 

da  — t 


and 


(3.12) 


— [wd] 

3t  1 V 


,2 

<>  , a,  , a,  a a a 

a — — (w  ) + c[w  ] = -c  w + f - cy  - cy  - eg 
3x2t  t t ttt  t 7 


w3(0,t)  = 0,  wa(l,t)  = g'(t)  , and 
t t 6 


w^(x,0)  = f (x,0)  - c (x,0)g^ (x)  + ag^(x) 


and  y satisfies: 
Jt 

(3. 13) 


Jt  [\]  ~ a 7T  [<]  " ° 

dX 


yf(0,t)  = g’(t),  ya (1 , t)  = 0,  and 
t 5 t 

ya (x,0)  = 

= 0 . 

By 

applying 

the  maximum  principle  one  easily  shows 

that  for 

all 

(x,t)  e 

(0,1)  x 

( 0 , T 1 

(3.14) 

ya  > 0 

(using  ( H 2 ) ) 

(3.15) 

ya  1 (1 

- x) [g  (t)  - g5<0)l  (using  (H2)) 

(3. 16) 

a 

w <0 

t - 

(using  (111),  (113),  (H4)  , (H5),  (H6)  , 

(H7) , and 

(3. 15) ) 

and 

(3.17) 

f^>  o 

(using  (3.14)). 
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- 


K(«t*t  ' *'  U'V  (HU)  ,m.l  (II'*)),  t fol town  t mmotl  t .1 1 n 1 y 1 1 om  t ho  nto.in  v»iluo  thooioin. 

Wi>  inm.iik  th.it  whon  t * 0 , o s 0,  mnt  >i  f ■ 0,  wo  .no  oxaot  ly  i i\  t hr  'ittu.xtton 

.m.il  y ■ i'il  in  | 'I.  wIkmo  tlio  tomtit  ot  Throi  out  1 fot  tin'  .ilvivu  !i|>ori.tl  mio  w.i  n tti-.t 
. 't't  .i  i nr. I . 

in  tin-  ili-i  t wit  i.'n  ot  tin'  i'iti'1  i"it  i rn.it  o wo  nti.ill  ini]ttttn  ,»  |'t  i v><  i nppm  .tint 
l.'Wi't  I'i'un.t-.  t,.i  tin-  mo  tut  ton  .1.  Wo  now  ptovoi 


-It 


Lemma  4 : Assume  the  conditions  of  Theorem  1 are  satisfied.  Then  0 < a4  < a <_  a* 


C t*  g'  (T)  1 2-V 

a*  = max]  1 , q /(/  5 dr  - g'(0))  \ 

1 L 8 0 Alt*'"-  t)  7 ' 


aA  = mi 


n{i'  KAlC  + g5(t*)  - Vt#,H4} 


(and  K is  a constant  defined  in  lemma  1). 
b 

Proof:  in  the  proof  of  Theorem  1 we  showed  that 

g8  = Q(a)  = -aya(0,t*)  - awa(Oft*)  - ag^(O) 

and  that  w (x,t)  < 0,  (x,t)  e I . Since  wa(0,t)  =0  it  follows  that  -awa(0,t*)  > 0 

~ T x 

and  so  g > - ay  (0,t*)  - ag’ (0).  Now  from  (3.18)  and  (3.19)  we  have 
8 — x 7 


and  so 


a ffc*  g (t) 

-ay  (0,  t* ) > a / 2 dt 

0 /na(t*  - t) 


If  a > 1 then 


r~  gMT> 

9 > A / ... . ?- dt  - ag'  (0) 

o *Ai(t*  - t) 


— ?5(-t)-  dT  - g'  (0)1 

' o A ( t * - tT  / > 


and  so 


( r ..  t*  q (!) 

a £ max  j 1 , g / ( J 

1 L 8 o A(t*  - 


-■  dx  - g^(0))j‘|  = a* 


To  get  a lower  bound  we  observe  that  from  (3.4)  and  (3.5)  we  haw 

gu  < 4 II  IF  - cza  - za)  ( • , t*)  II  + a[g  ( t* ) - q (t*)  ) < 4 Kr  + a [gr  (t*)  - q ( t * ) 

b — j ' t U 5 6 ~ 3b  5 b 

for  0 < a < 1/4  (by  lemma  1) . Hence  for  0 < a < 1/4 


and  so 


"ail!  '? 


“ hs/lK*  " V 


t*)  - Cl  (t*) 

6 


D4)  = *.  • 


. DESCRIPTION  OF  THE  APPROXIMATE  PROBI.KM  AND  ERROR  ESTIMATES  FOR  PROBLEM  (P) . 


It  is  convenient  to  descirbe  the  approximate  problem  in  somewhat  abstract  terms 
since  part  of  the  numerical  procedure  involves  the  approxim.it  ion  of  the  forward 
problem  (l.S)— (1.7)  for  which  many  numerical  schemes  are  available. 

For  a given  lot  lia  denote  an  approximate  solution  of  (l.fi)-(1.7)  corresponding 
to  meshes  of  width  h and  At;  in  the  x and  t directions  respectively  and 
satisfying  an  error  estimate  of  the  form 

(4.1)  lu^O.t*)  - lfl'(0,t*)|  < C (b  , b*  )<£  (li.  At) 

x x 1 — * 


for  all  h c |b.,b*]  where  ^>(h,At)  satisfies  limsf(h,At)  = 0 (for  some  schemes 

h Hi 
At  HI 

this  type  of  estimate  may  only  hold  when  some  relation  exists  between  the  size  of  h 

2 

and  At,  e.g.  At  Ch  ) . 

The  approximate  problem  is  then  to  determine  a constant  « such  that 
(4.2)  C>(u)  s-au”(0,t*)  = gg  . 

Since  this  is  just  a nonlinear  equation  in  a,  wo  can  solve  it  by  any  standard 

it 

iterative  method  where  of  course,  in  order  to  evaluate  Q(«  ) we  must  find  U , 

the  approximate  solution  to  the  forward  problem  (1.5)— (1.7)  (with  a replaced  by  <r^) . 

Wo  remark  that  for  finite  difference  scheme*;  tho  derivative'  U (0,t*)  may  be 

x 

$ 

replaced  by  a difference  quotient  such  as  |U  (h,t*)  - U ( 0 , t * ) ) /li , and  Theorem  2 
(to  follow)  easily  amended  to  account  for  this  change. 

For  approximation  schemes  satisfying  (4.1),  we  will  first  show  the  existence  ol 
a solution  to  the  approximate  problem  and  derive  an  abstract  erroi  estimate.  We  then 

i 

give  specific  error  estimates  for  two  examples  where  an  estimate  of  the  type  (4.1) 
can  be  easily  derived  from  results  appearing  in  the  literature. 

Theorem  2:  Suppose  tho  hypotheses  of  Theorem  1 are  satisfied  and  that  ( a , u'  ) is 
the  solution  of  problem  (1.5)-(1.8).  For  a given  let  U‘*  denote  an  approximate 
solution  of  (1.5)  — (1.7)  satisfying  condition  (4.1).  Let  0 < b^  *-  a4  «-  a*  < b* 
where  a#,a*  are  the  bounds  on  a derived  in  lemma  4.  Then  for  h and  At 

\ 
\ 

* 
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sufficiently  small  there  exists  an  a r |b<fb*)  such  that  the  pair  (a,lia)  satisfies 
(4.2)  and  for  all  such  solutions  we  have 

<4-  -»>  |a  - a|  •_  c(b*.b*)v*  (h,  \t  ) . 

Proof ; By  definition, 

0(b)  - -b  Ul*(0,t‘)  -b  ul’*(0,t*)  t b . | iil  « (0, t *)  - l'b*(0,t*)l 
x * x * x x 


0(b4)  ' C (b^ , b*  (h,  t)  (usinq  (4.1))  . 


Now 


Q(b.)  ■ 0(b  ) -Ola)  ♦ q < q - K(a*)(a  - b.) 

" 8 H " 

(usinq  ( 1.20)  and  the  fact  that  a^  «-  ,i).  Hence 

01b,)  ^ q(,  ~ lMa*)(a  - t>4)  - 0 (b^ , b*)^  (h , At ) 1 *■  for  h,At  sufficiently  small 

since  C(h,At)  *■  0 as  h,At  * 0.  Similarly, 

Q (b* ) > 00’*)  - 0(b#,b*)vMh,At)  q + K (b* ) (b*  - a)  - C(b#,b*)<  (h.At) 

N q for  h,At  sufficiently  Small. 

By  the  continuity  of  0(a),  there  exists  an  « < lb4,b*l  such  that  0(a)  q, 
for  all  h,At  sufficiently  small. 

To  derive  the  error  estimate  we  choose  b a in  (<..’)  to  obtain 


I*  - “I  - K (b *T  ^(a)  " «(o)l  1 K • K K T(i*)  l^(a)  - ^ 


l 


- 0(a) 


b* 


- zr^iT  U (0,  t * ) - u (0,t«)  v C(b  »b*)V  ()i,At) 

Mb*)  1 x x 1 — Mb*)  * 

(usinq  (4.1)  and  the  fact  that  a < (b#,b*|). 

A;,  an  application  of  Theorem  we  shall  consider  the  Crank-Nicholson-Oalerkin 

approximation  to  (1.5)— (1.7),  which  may  be  described  as  follows:  bet  h » IN,  where 

N is  a positive  integer  and  sot 

sf  - Iv.  C°|0,  U : V « P (I.)  1 1 N)  , 

h k - 1 i , 

where  l • (li  - l|h,ih)  and  P.  . (1)  is  the  set  ot  nil  polynomials  of  deqree  v k-1 
1 * “ • 

defined  on  1.  Define  sj4  (v  < S,  : v(0)  Vil)  Ol.  Now  let  At  - T where 

h h 

is  a jxisitive  Integer  and  set  t ■ nAt,  n 0,l,...,J. 


For  a function  F defined  at  times  t we  denote  by  F the  function  Fit  ).  We 

n n n 


shall  also  use 


F , “ (F  ♦ F )/2 

n+1/2  ntl  n 


(Note  F * F (t  . . 1)  and 

n+1/2  n+1/2 


3 F - (F  - F ) /At  . 

t n«-l/2  n+1  n 


Our  approximate  problem  is  then  described  as  follows: 

Determine  a constant  a and  a sequence  < lirt  > where  ll"  - za  + r , Zn  f S** 

I n / n=0  n n n n 1 

satisfies 

l4'4’  ('tCt/rvl  ' •"£  Cl/2’  ir*  vl  * es-.m.i|*oz“tl/J.v>  - If„,1/2.V) 

ck 

tor  all  V f S,  , n > 0 and 
h — 

, . , . ,3  , a ..9  . _a  «k 


J<3x  [Z0  " G| ' * <c<-,0)IZ0  - G] , V)  » 0 for  all  V r Sj 


whete  r » (1  - x)g_(t)  + xq  (t)  , 

5 6 

F «»  f - cr  - r * and 
G - q.;(x)  - r (x, 0)  , 

and  - all  (0,t*)  - q . (z  is  the  Crank-Nicolson-Galerkin  approximation  to  s'), 
x 8 n 

For  the  case  when  c - c(x),  Wheeler  [15|  has  proved  that  if  z (defined  by 
k rt'  3 7 k 3 7 

( i.  1)  1 belonqs  to  1.  (W  ’ ) , — belongs  to  L (H  ) and  — - belongs  to  t.,(l.,), 

3t  2 3t3  - 2 

then  there  are  constant:.  0 and  t ^ 0 such  that  for  all  0 < At  < t „ 

0 0 


U.M  !*"(•. t„l  -z“l  iclhk(|UHL  Hllfll  ,„•>,>*  (AO 


•’  tl  3 7.  I 


3t3"  W 


, n "*  0 . 


Although  the  constant  C will  depend  on  a it  is  «.  c(b#,b*)  for  n r (b^ ,b*| . 
Remarks:  In  [14]  Wnhlbin  showed  that  an  estimate  similar  to  (4.6)  could  be  derived 

for  moie  easily  computable  choices  of  Z(1  than  (4.5)  and  also  for  subspaces  with 
higher  interelement  continuity.  Wo  further  note  that  it  is  possible  it  extend  the 
results  in  [15|  to  the  case  c - c(x,t)  for  the  problem  wo  are  considering.  (A  pioot 
of  this  fact  is  included  in  the  appendix  as  lemma  6) . 
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To  get  an  estimate  of  the  type  (4.1)  we  use  an  inverse  condition  satisfied  by 
k 

elements  V * s,  , i.e.# 

n 

(4.7)  llv.ll,,  -f  I|V|IL  ' 

GO  00 

oi  k 

Letting  V be  the  interpolate  of  z (*,t  ) in  S,  , we  get 
I n n 

Ilf  [ua(.,t  ) - ua]||T  - Ilf-  [zV.t  ) - Za)|| 

1 3x  n n ' L 9x  n n L 

00  00 

l*°<-.tn)  - *,1||L  *fllvi-Z«||L 


ills  - v«L  ♦!»*,  - - z“„i'l 

00  00  oo 

i»k‘1||*a<-.t  )||  tflliV.t  > - z“|lt 

W 00 

k 

(using  the  standard  approximation  properties  of  subspace  S ). 

h 

k-1  2 

Applying  (4.6),  we  will  then  get  the  desired  estimate  with  'f  (h,At)  = h + (At)  /h 
provided  that  the  solution  za  has  the  indicated  regularity  for  all  a e [b^,b*] . 

Using  Theorem  5.3  in  Lions  18),  this  will  be  the  case  if  the  following  compatibility 


conditions  are  satisfied. 


G(x)  =0,  x = 0,1  , 


(4.10) 


aG"(x)  + c(x)G(x)  + F (x, 0)  =0,  x = 0,1  , 


a [aG" (x)  + c(x)G(x)  + F(x,0)]" 
+ c(x)[aG"(x)  + c(x)G(x)  + F (x, 0) ] 


+ Ffc(x,0)  = 0,  x = 0,1  . 

Since  we  require  that  conditions  (4.9)  and  (4.10)  hold  for  all  a e [b^.b*] 
these  two  conditions  actually  lead  to  the  more  restrictive  conditions  that 


(4.11) 


(4.12) 


(4.13) 


(4.14) 


G" (x)  =0,  x = 0,1  , 


c(x)G(x)  + F (x,0)  =0,  x = 0,1 
G(lV>  (x)  =0,  x = 0,1 
(c(x)G(x)  + F (x,0) ] " =0,  x = 0,1 


and 


(4.15)  F (x,0)  =0,  x = 0,1  . 

(For  the  case  c = c(x,t)  the  compatibility  conditions  become  more  complicated.  See 
Theorem  6.2  of  (8)). 

Under  these  conditions,  we  get  the  error  estimate 

| a - a | <_  C(b^,b*)  [hk_1  + (At)2/hl  (k  < 4)  . 

By  making  use  of  the  results  of  11)  and  observing  that  condition  (4.1)  only 

involves  an  estimate  at  time  t = t*  > 0,  the  compatibility  conditions  stated 

previously  can  be  weakened  for  a special  case  of  the  problem  we  are  considering. 

Suppose  that  g_  and  g,  are  constants,  c = c(x),  and  f = f (x) . Let  o'  (x)  be 
5 o 

the  solution  of  the  steady  state  problem 

.2 

(4.16)  -a  — — 0 +c0  = f,  0<x<l 

dx" 

0a(O)  = g , 0a(l)  = g . 

3 b 

a a a 

Then  v = u - 0 satisfies: 

(4.17)  va  - av*  + cv°  =0,  0 < x < 1,  t > 0 

t XX 

va(0,t)  = 0,  vn(l,t)  =0,  t > 0 
va(x,0)  = g?(x)  - 0a(x) 

d 

Now  denote  by  •6*  e the  standard  Galerkin  approximation  to  the  steady  state 

solution  0a  and  by  {va}  the  Crank-Nicolson-Galerkin  approximation  to  v". 

n 

It  is  shown  in  [1]  that 


(4.18) 


|v  (* ,tn)  ‘ n I l — c (tn> lh  + <At)  nlg7-  e II o 


provided  At  <_  Ch  . Since  we  have  from  [16]  that 

k||  .Oi 


(4.19) 


|0“-^|L  < chK||ea||  kf-  , 


it  follows  that  for  U 


,a 

n 


Va  + *6? 
n 


•! 


I 


1 
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C(h,.b*)|hk  * (At)’’] 


|uV.t*)  - Ua(-,t»)  lk  i 

OB 

k—  1 2 

and  using  the  inverse  condition  as  before  that  (4.1)  holds  with  <f  (h, At)  li  + (At)  /h 

without  requiring  the  corresponding  compatibility  conditions: 

(4.20)  I g " ( x ) - 0 " ( x ) ) - 0,  x * 0,1 

and 

(4.21)  c'(x)(q'(x)  - 0*(x)l  =0,  x “ 0,1  . 

2 

Choosinq  S the  space  of  continuous  piecewise  cubics  (i.e.  k 4)  and  At  = Ch*-, 
we  would  then  qet 

(4.22)  | a - a | - c(b(,b*)h  ' . 

We  note  tliat  results  analogous  to  (4.18)  are  obtained  in  (11  for  a wide  class  of 
single  step)  Galerkin  approx  un.it  ions , and  that  error  estimates  corresponding  to  use  of 
these  schemes  are  obtained  in  a similar  fashion. 
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A A 

Now  from  (5.2)  we  have  -aAz  ■ f - ixj  - cr  and  so 


a||z'*|i  , - ||- A 1 ( f - cw  - cz'‘|||  < ||- A 

9 ’ U)‘  * 


!I-a-'|I  ,,  . II f - « - c*|L 

« I 4 

— 4 0 


by  the  Sobolev  lemma  — lx*)  v c z ^ . and  so 

3n  ” ..2,4 

w 


a (x*;  I < c||f  - CU)  - «a||L 


It  then  follows  that  l im  Q(a)  0 since  ||  f - cu>  - e; 


| ► 0 as  a * 0 by  lemma  7.  . 

^4 


3 

multiplying  equation  (5.2)  by  [z' ] and  integrating  by  parts,  wo  have 
3ai|zAVz,H”  + ||/c  [z'}  (1 1 - (f  - cui,[za]  ) 


< f - Cii» 


lL  IIUa]\  -llf  - mK  ll*allL 

4 4/3  4 •! 


~ a 

Hence  | ) z ^ 1 1 ^ — ||  f - cui||  (using  A3).  This  implies  that  |a  -r— — (x*)|  is  bouiided 

L4_C1  L4 

tor  all  a and  thus  from  (A2)  and  (5.3)  it  follows  that  lim  Q(a)  = +*“.  The  existence 

a *«* 

of  a positive  solution  to  i^(a)  g*  for  [xisitive  g*  follows  from  the  continuity  of 


To  prove  continuous  dependence,  we  differentiate  Q (a)  to  obtain 

3 3 A 3 , 

Q’  (a)  = -r-  [a  ( X* ) 1 + ~ (x*)  . 

3a  3n  3n 

Now  it  is  not  difficult  to  show  that 

3 3za  3 3za  3za  3 3za  a 

— (a  lx*)  1 = a — - (x*)  + - — (x*)  = -r—  la  ^ — + 2 ' 

3a  3n  3a  3n  3n  3n  3a 


3z 

where  — — satisfies  the  boundary  value  problem 

i’ll 


- Aza  i, 


0 on  P . 


Setting  p*'  = a ——  + z'',  we  have  that  O’ (a)  = ( x * ) + ~ (x*)  where  p satisfies 


.a  a .a  , a , a , 

aAp  ♦ cp  = aAz  - aAz  *■  cz  in  s; 

pa  = 0 on  r 


r 


Applying  the  maximum  principle  it  follows  from  (Al)  and  (5..))  that  «■  0 and  hence 


- 

ap 

, r4—  < x * ) ' l and  ■ <> 

3n 


from  (5.5)  that  pa  - 0 in  (1.  Since  F'1  = 0 on  r 
O' (a)  2 y~  (x*).  (5.1)  now  follows  directly  from  the  mean  value  theoiem. 

As  in  tlie  parabo lie  problem  we  shall  require  a priori  uppei  and  lowet  bounds  foi 
the  solution  a in  order  to  derive  the  erroi  estimate.  We  now  prove: 
lemma  5:  Assume  the  conditions  of  Theorem  t are  satisfied.  Then  0 » a 


a < a*  for 


some  constants  a , a*  depending  only  on  the  data  ot  Problem  (F.) 
P t oo f : In  the  proof  of  Theorem  3 we  showed  that 

g*  = 0(a)  m a (x*)  + a ~~  lx*) 

3n  3n 


and  that  za  v 0 in  i).  Since 
3u> 


a 


(x*).  hence  by  (A2)  a < g*/ 


0 van  P,  it  follows  that  (x*)  > 0 and  so 

3n 

3u> 


” 3n  ~ ” 3n 

To  vjet  a lower  bound  we  observe  that  from  (5.3)  and  (5.4) 

g*  <c||f  - cw  - c*a||L  + a (x*)  < cal  + a (x*) 

4 

for  0 < a < (bv  lemma  2).  Hence  for  0 a (*  t defined  in  lemma  2) 

0 0 0 

1/8 

< Ca  and  so 


a N min{e_#  [q*/cr  1 « a*  • 

Remarks : We  now  make  two  comments  about  the  requirement  in  (A3)  that  six'  ''  c ^ N 0. 
As  evident  from  ttie  proof  of  Theorem  3,  the  reason  for  this  restriction  was  tv'  prow 
existence  of  vi  solution  for  all  positive  vi* . Wo  note  however  that  it  c only 
svitisfies  c(x)  2 then  it  is  still  true  that  O'  (b)  N 0 for  b '*  0.  Hence  it 

we  also  know  that  - g*  and  0(ao)  ~ q*  witti  0 < a *-  a f then  for  .my 

.I*  v q*  2-  g*.  there  will  exist  an  a^  ^ a ' a,  such  that  0la)  g*  • Since  we  would 

thi'ii  have  both  the  existence  of  a solution  and  a priori  bounds  on  a without  the  use 
of  lemma  2,  the  requirement  that  c(x)  2 '•'j  N 11  could  be  weakened. 

The  second  remark  lias  to  do  with  the  case  cixl  _ 0.  In  this  case  :■  a 

wtiere  ?.*  satisfies 


- Az 


f in  U 
0 on  F . 
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to.  DESCRIPTION  OF  THE  Al’PKOX  I MATE  PROBLEM  AND  ERKOR  ESTIMATES  1 OK  PKPHLEN  (E)  . 

As  in  Section  -1,  it  is  convenient  to  describe  the  approximate  problem  in  abstract 

terms  so  that  we  may  easily  take  into  account  the  many  numerical  schemes  available  for 

solving  the  forward  problem  ( 1 . 9) - ( l . 10 ) . 

a 

For  a .jiven,  let  U denote  an  approximate  solution  to  (1.9)- (1.10)  satisfying 
an  error  estimate  of  the  form 

(b.l)  |—  (x*)  - (x* ) | < C (b  ,b*)y(h) 

an  an  — * 

tot  all  (b<(b*]  where  y (h)  satisfies  lim  y (h)  = 0.  (In  the  applications,  h 

h *-0 

miglit  denote  the  size  of  a finite  element  mesh) . 

The  approximate  problem  is  then  to  determine  a constant  a such  that 

(to.  2)  Q(a)  = a (x*)  = g*  . 

Once  again  this  is  just  a nonlinear  equation  in  a,  easily  solvable  by  standard 
iterative  techniques. 

Using  essentially  the  same  proof  as  in  Theorem  2,  we  get: 
la.iem  -1  suppose  the  hypotheses  of  Theorem  3 are  satisfied  and  that  (a,uJ)  is  the 
.,'luti  >n  of  Problem  (E)  . For  a given,  let  Ua  denote  an  approximate  solution  of 
(1.9)  (1.10)  satisfying  condition  (6.1).  Let  0 < b^  < at  < a*  < b*  where  a4,a*  are 

. rounds  on  a given  in  lemma  5.  Tlien  for  h sufficiently  small  there  exists  an 
■ ' .b  ,b*  mdi  that  the  pair  (a, ua ) satisfies  (6.2)  and  for  all  such  solutions 
we  have 

(t  . t)  |a  - a|  C(bA,b*)y (h)  . 

in  order  to  apply  Theorem  4,  we  need  only  choose  the  approximation  scheme  and 

determine  the  function  y (h)  of  (6.1)  Estimates  of  this  type  for  various  finite 

element  methods  follow  easily  from  the  work  of  Nitsche  (9],  Goldstein  and  Scott  [6J, 

Wahlbin  113]  and  others  described  in  the  references  of  those  papers,  provided  we  have 

the  inverse  condition  II W II  < C-  ||v||  for  all  V e finite  element  subspace 

L — h !■* 

OP  CO 


beinq  used. 


7.  APPENDIX. 


Proof  of  lemma  1;  For  any  — > e > 0 define 


X (x)  = 


x/e 

1 


0 <_  X <_  E 
e < x < 1 - e 


Then  one  easily  verifies  that 

(7.1) 

(7.2) 

(7.3) 


(1  - x)/e,  1 - e <_  x <_  1 . 

o i xc  1 1 

c 


1 - X 


</n 

0 - * 3 


and 


l^xe||0 


a a a a 

Write  z ~ zi  + z2'  where  satisfies 


(7.4) 


9 a 92  a 

9t  Z1  ' 3 


Dx 


d 3 G O 0 

2 Z1  + CZ1  = X fzt  + CZ  1 


and  z^  satisfies 


(7.  5) 


za(0,t)  = 0,  za(l,t)  = 0,  za(x,0)  = G (x) 


9 a 92  a , a e. 

91  Z2  " 3 TIT  Z2  + CZ2  = FU  ' X 1 

3x 


z3(0,t)  = 0,  z3(l,t)  = 0,  z3(x,0)  = 0 . 


From  (7.5)  it  follows  that 


,9  a a.  , .9  a 9 a.  , a a. 

(9l  W + 3(97  V 97  Z2)  + 

= (Ffl  - xE],z3)  < ||  F ( • , t ) ||  ||za(-,t) 


(using  (7.2)) 


Hence 
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Applying  Gronwall's  inequality,  we  get 


Ho 


i # I / o'"8  ||  F ( ‘ Is)  ||f  d>) 


1/2 


Now  using  (7.4),  it  follows  that 
3 , a c 0.  a i 0 


, 3 . a c u a to  ,3  a 3 a,  . . a t u a c u 

(^[*1-XM,*1*xO+  a(-  *1#  - zL)  + (c|Zj  " X * " X * ) 


f 0, 


C 0. 


,3  a 3 , e 0 |i  3 a . ■ ■ ’ 

= c lx  2 n < aIUr  z,  (-.t) 


3x  1 3x 


3x  1 


0 4 11 3x 


Hence 


JL 


e„0. 


ii  »IZI  - x'*'n-a)||5  if  ||£  ixEo°e,t)i 


lo 


and  so 


, a i 0 . . 1 1 2 

[Zl  - X 2 I (-ft)  ll0 


[1  - X'JG(X)  L t 


|/  ll£  lX,2°(.,S),||f)dS 


Now 


& z> 'll oi II >S  ll-u‘-.  = » *X‘II0 


c 3 0 


3 e, 


- z°(-'s)  II 0+  J1  ||z°(-,s)  ||  L 


(using  (7.1)  and  (7.3)).  Hence 


lz^  - z°l  ( • ft)  ||~  < 


V?  i|u||l 


1 2 


t 


+ a / 
0 


3 °,  x||2. 

z C,s)  ||0ds 


+ a / 
0 


7 2 


ds 


and  so 


[z*  - Xlz°l  *fc)  ll0  i V“T  llGllL  + ^ i f 

0°  0 


t 


1/2 


+ 2 | /^  {/  ||  z°  ( * ,s)  llj'  dsl  1/2 


Applying  the  triangt  inequality,  we  get 


[za  - z°I  (*  ,t)  ||0  1 


[z^  - xrz°H‘»t)  ||0  + II Z2  < * ' *")  ^1 0 


+ ||  tXC  - 11  z° ( - ,t)  || Q i kJ/T  + \/f 


wtier  o 


Ki  l/"r  llull,.  * */  °l  ’ll  k ( - , s>  || * dsi  + ||z°(-, 


• nut 


t 


K 


1 


</’  II  S '°<-.»lloJs|1/i 


t n > 1/2 

21/  ||  7.  (•,a)||“  ds} 

0 


Choosing  i /a,  wr  have  for  0 < a < 1/4  that 
(7 . d) 

when1 


a 0 . ||  t 1/4 

z - z ] (*,t)  ||  0 v K^a 


Ki  * K1  + ^ + K3  • 

To  obtain  an  estimate  for  |||/.'('  - z^]  ( • ,t)  ||()  wo  observe  that  since  G*  c“lO,  1 
with  GtO)  G(J)  0,  /.'*  satisfies  the  initial  boundary  value  problem 


2 a 
9t  :'t 


2 ‘ a a 

a z.  • cz  - F.  - c z 


9x 


t t 


and  z sat  is  t l e 

t 


z‘  (0,t)  0,  7.‘  (1  , t)  0 


(x,0)  aG"  (x)  + F ( x , 0 ) - c (x,0)G (x)  , 


9 0 0 t 

..  7.  * cz  • F.  - c z 

3t  t-  t t t 


z (x,0)  K(x,0)  - c (x 1 0 ) G ix ) . 


Wo  now  apply  the  previous  argument  to  those  problems,  i.e.  writ* 


t 

9 <i 
li  t ;’t  ,1 


a a a a , , . 

7.  . t 7.  whoio  7.  , satisfies 
1 1 t t 1 


(7.7) 


9 a , a i . 9 0 

- ,1  v»  ♦ c?.  = V I — 7.  ♦ 

„ 2 tl  tl  x 1 9t  t 

9x 


Z*  , (0,t)  0,  Z . (1  ,t  ) 0 


1 1 


1 1 


z (x,0)  aG"(x)  * F ( x , 0 ) - c(x,0)G(x) 
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and  2t.,  satisfies 

> 9 a 

(7-8)  9l2t2 


92  a 


a 


a — * t czt2  - U - X 1 lPt  " ct*  ] - ctl«  - 2 


z'\,  (0,t)  0,  z*  (l,t)  = 0,  zf  (x,0)  - 0 . 

t2  t2  t2 


Proceed i ru?  as  before,  wo  obtain  from  (7.8)  that 


{/  et  " \Py  lhFt  " ctz°  (-'s)ll^  ds} 


S2  i/  el  "lie.  (.,s)||;  III/.'1-  z°]  (*,S)  ||  ‘ ds) 


/Via in  following  out  previous  argument,  we  get  using  (7.7)  that 

Hl\r-  X*  z^l  ( • , t ) ||  Q 1 a ||G"  ||  0 + ||F(-,0)  - c ( • , 0 ) G | 


/a  { / 


/ II  ;j~  zt(‘ ,s)  II 0 ds)  ' 2\[~  II ( * . s)  II L ds} 

0 0 <•> 


Applying  the  triangle  inequality, 


Iztl  “ xL*t1(*,t)  Hq  + Ilzt2(,'t)  Ho  * Hlx1  - l)z^(-.t) 


Choosing  t = »a,  inserting  inequality  (7.0),  and  collecting  terms  we  will  have  for 

some  constant  K*”  that  for  till  0 < a 1/4 
5 

1 1 , a 0 _ . -II  t 1/4 

II  (zt  - zt]  (‘ ,t]  II 0 s K5a 

Inequality  (2.3)  now  follows  from  the  definition  of  z°  and  the  triangle 
inequal ity . 


1'roof  of  lemma  2:  For  1 s 0 define 


XC(x) 


1 , p ( x , r ) s t 


p (x , r ) , p ( x , r ) < 


where  p (X,  n denotes  the  distance  of  a point  x e il  to  the  lioundary  T.  Since  C 

is  assumed  sufficiently  smooth  there  will  exist  an  1 1 > 0 such  th.it  for  all 

0 

r 1 ,4 

0<ev_r  , x ( W and  satisfies 


(7.9) 

0 < xe  < 1 

(7.10) 

111  - XE|Il  £Ce 
L4 

and 

(7.11) 

I|VXC||L  £ Ce" 

1/4 


-3/4 


^ a 3 a . a . _ . 

Now  write  z = Z1  + Z2  w^ere  satisfies 


(7.12) 


3 c 

(-aA  + c)z  = (f  - cw)x  in  ft 


z^  = 0 on  T 


and  z^  satisfies 


(7.13) 


(-aA  + c)z^  = (f  - cu) (1  - xC)  in  ft 
z^  = 0 on  r . 


From  (7.13)  it  follows  that 


-a^z*,^3)  + (cz*,lz*)3)  = (If  - coj 111  - x^U*]3) 
and  after  integration  by  parts  that 

3a||2*Vz*||^  . ||c1'V||4  < ||  f - »||  ||1  - x'  II.  ||l«*)3|| 

4 °°  4 4/3 


< C||f  - cw||T  e1/4  ||  zj  ||^ 


(by  7 .10) ) . Hence 
(7.14) 


(since  c(x)  ^ c > 0) . 
Setting 


I Z2  II L - CHf  " C<1)HL  cl/4 

4 CO 


.£  rf  . E 

'I'  = [-  - ui]  x , 


it  follows  from  (7.12)  that 


-a(A[z*  - *E],Iz*  - *V)  + (c[z*  - *E]f[zJ  - i/>E]3)  = a(Ai).e,fz^  - *E]3) 
and  after  integration  by  parts  that 
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I U --  1 11  1 


i«t  - *'iv»'n0  hi,;  » i n0 


,J"ih<  - ♦•iv.’i*  • ...in/;  - v'ivi i n- 


Monet* 


III/' 

4.  ■ 11  1 


/ 1V*‘  II ;;  v II /;  - <•’  II;’  IIV  ||; 


■ Sllvi  SO 


( 7 . I *>  ) 


11/  <■'  II,  • s//  II  s 


Now 


(7.  UO 


Ill'll,  • IK'VI1.  - ...III,  • III'  ' ...  I ' 


Vl’  ...III,  ' /II, - ...II,  ."  '' 

•1  'v 


(by  (7.11  M.  M'liw  1 1*«  * t t i .»n«}  1 « * im**|n.il  »l  y wo 


1 1 - v*<i>  - cr 


V-  ,1  I)  l - U - \ l 


• IU-*  - /' 


i "i 


/i. 


, l 4 1/2  -44  l 

C 1 1 * .1  » i .1  I 


looml' l ni  n>)  out  previous  tesults  l ’.14)  i'Iu'oh  uw 


Wo  luvo  t.'l 


0 *•  .1  *■  i tls.st 

0 


I - o«>  - O 


1 


for  some  eonsitnist  0 i no  only  on  t,  >■ , ut  .uv.l  s'. 

I.emiss.i  us  t'Upixrso  that  ••  is  t In-  sol\»t  ion  ol  rroblom  (<.*'  .nut  l.-'  < >'«'  solution 


(4 . 4)  - (4 . . It  o(x,t)  0 e U.O  < 0 tos  ix.t'  « 1 1)(  -uui  .•  ss 


.sn  ti  i . ■ i out  1 y smoot  Is , t lion 
(7.1 7) 


St  ( • , t I V S’ 

n is 1 1 


.>  , 


k ,.v,  ■ ,v i 

l (w  i l ,01  i 


(At  > 

\ 

• „ , * 

ll«ll 

\ * 

1 ,0*  / 

M 

t.  . (i.  / 

for  some  constant  C independent  of  h and  At. 


Proof : Following  the  proof  in  Wheeler  [15]  we  define  for  every  t f 10, TJ,  z(',t)  e S * 


(7.18) 


3 - 3V  ~ 

a <]T-  (z  ~ z]  ,r~)  + (c[z  - z]  ,V)  =0,  V e S,  . 


3x  - 3xr  • ” ’ - “h 

Averaging  the  weak  form  of  differential  equation  (3.3)  at  times  t = t and  t = t 
yields 

3 3 


(7.19) 


(3  z 
= (F 


t*n*l/2,V’  * Vl/2'ta  V|  + l,clln.l/2'V) 

“k 


„*1/2’V1  * <VV>-  VtSh’ 


n.  1 


where 


1 n+1  3 3 

' n ^tZn+l/2  ^Zt^  n+1/2  “ 2At  ^ ~3  (t)  (tn+l  ~ t)  (tn  ~ T)dT  ‘ 

t 9 1 

n 

Now  setting  n = z - z and  £ = Z - z , we  have  for  all  V e s?4  that 

n n n h 

(9tW2'V)  + a(^T  Cntl/2'4  V)  + (C(-'[n  + I)At)Cn+l/2'  V)  = "(£n  + 3tVl/2'V) 

+ (fCZ)n+l/2  " c(-'[n  + I]At,Vl/2'V) 

(using  (7.18),  (7.19),  and  (4.4)).  Setting  V=  3t^n+1/2>  we  get 

<7-20'  “Vn.l^O  * m - ll£  «Jo>  . l|cV2(',[n  . i]At,Vl||2 


-ll=1/2c,ln  . iut.cjl2  <i||en  . atW2l|2  .||atc 

* 7 11  ,c£,„.1/2  - * JI4t’V1/2llo  • 


n+1/2  11  0 


Since  c < 0, 


|c1/2C.ln  . iMt)£n41||2  > ||c1/2(..[„  ♦ flAt)£ntl||2 


Making  this  replacement  and  summing  inequality  (7.20)  from  n = 0,...,J*  - 1,  we  get 

ll^50llo!  ♦ll<=1/2<-.w*.  jHet  .llo -llc1/2i-.jAt)E0llo 


.ii; 


J 


* 

J -1 


12  _L  {||cn  + 9tr'n+l/2  ^ 0 +Hlczln+l/2  " c(''[n  + ? 1 At)  Vl/2  Hq} 


n=0 
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Since  = 0,  we  have 


\\h  < . 


IJ,  - .Jo  i l <11%  * VW2H0  * J|C2l„.l/2  ' Cl’'ln  * 2lati;n.l/2llo< 

J n 0 


J -1 


Now 


and  so 


n+1/2 


[c  4 c 1 

ntl  n 

- ■ 

z , 4 z 
n*-l  n 

4. 

c . - c 

n+1  n 

r-  * “ 

z . - z 
n + 1 n 

2 

2 

. 

» 

j 

1 ; 

NJ 

1 ! 

2 

_ _ 

l'lw,nU/2  " c('.(n  ' 2 <Atl/'n»  1/2^0 


< llVilla  - c ( • , tn  . i|4t)  II  t llv1/2ll0  * ll-aiiT~al 


z , - z 
1 n+1  n 


-(  4 ^Ctt^1»(L®)  ^*n+ 1/2^0  * 4 ^'tZnU/2^0 


Estimatinq  terms  as  in  ( 1 S 1 we  qet 


(7.21) 


1st e -»5  if  { 

J 


(At) 


33z 


At 

2 

At 

2 


3t  3 
2 r 


W 


I— I 

’at1 


W 


tt"l,  (L  ) 


2HZIIl2(L2)  + (At)  'Ut^L2(L2)J 


t L (L  ) 

00  on 


0*2 

I— II 

'3t"L2(L2 


)}  — C < * » a ){ 


(At) 


a3 
a z 


at' 


2 


w 


♦II^cl,.  * K2,l2,  * «*•' 


WIVV  ♦»H«l2,lj, 


Iq  1 C(a^,a*)hk  ||  z ( * ,t)  Ilk 


(since  z = n 4 z) . Since  by  the  standard  Nitsche  arqument 
(7.22)  ||  n ( - ,t) 

we  will  qet  the  desired  result  by  following  the  arqument  in  [15]  provided 

+ lull  . 1 . 


ifjllj  (I  . 1 c(a  fa*)hk|||z  II 

3t  L2(L2)  h2(HK) 


k 

l2oi  ) 


Now  f rom  (7.18)  we  have 


"•  h v>  ♦ <cn'v>  • »•  v,lik,  • 
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Hence 


(7.23) 


To  estimate 


a(Fx  V £ v>  + (cVv)  = -<V'v)  . 


1 3n 


3t*'L.  (L  ) We  define  for  every  t e [0,T] , ij/(-,t)  e S* 


by 


(7.24) 


3 3 g 

a(3x  Izt  ~ ' 3^  V)  + (cfzt  “ «v)  = °'  V e sk 


h ' 


Then  from  (7.23)  and  (7.24)  we  have  for  all  V e Sk  that 

h 

a(^  Ut  -tl^V)  + <c[zt  - i|/],V)  = (-ctn,V) 
Choosing  V = z^_  - ^ one  easily  gets 

II  [zt  ~ !/-](-, t)||0  < C(a#<a*)  ||n(*,t)||  Q . 

Now  by  the  standard  Nitsche  argument 


(7.25) 


II  lzt  ~ *1  (*.t)  ||0  1 C(a^,a*)hk||zt(.>t)||  k . 
Hence  from  (7.22)  and  the  triangle  inequality,  we  get 

Hnt(,'t)  II 0 i C(a,,a*)hkr||zt(.,t)  ||k  + ||z(-,t)  ||k) 


and  so 


^ nt ^ l2 (l2>  — c ^a* ,a  ll|z 


f-  k + 11*11  k ) • 

l2(hk)  l2(hk) 


The  remainder  of  the  proof  of  the  lemma  is  identical  to  the  proof  in  (151. 
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A priori  error  estimates  are  derived  for  the  approximation  of  an  inverse 
probiem  in  which  it  is  desired  to  identify  an  unknown  constant  coefficient 
in  a one  space  dimensional  parabolic  or  two  space  dimensional  elliptic  partial 
differential  equation  whose  general  form  is  known. 
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